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Abstract. We study the direct and inverse spectral problems for semi- 



classical operators of the form S = So + Fi^V, where 5o = | ( — fi^ Am" + 



\x\'^ ) is the harmonic oscillator and V : R" — )■ R is a tempered smooth 



function. We show that the spectrum of 5* forms eigenvalue clusters as 
h tends to zero, and compute the first two associated "band invariants" . 
We derive several inverse spectral results for V, under various assump- 
tions. In particular we prove that, in two dimensions, generic analytic 
potentials that are even with respect to each variable are spectrally de- 
termined (up to a rotation). 
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1. Introduction 

Consider the semiclassical harmonic oscillator 

So = ^(^-h^Aun + - hnl^ 

acting on L^(R"), where, to simplify the notation throughout the paper, we 
have subtracted the ground state energy. The spectrum of Sq consists of the 
eigenvalues 

(1.1) Ej{h) = hj, i = 0,l,2,---, 
with multiplicity 

fn + j-l\ 
"^^■=1 n-1 )■ 
We will study perturbations of 5*0 of the form 

(1.2) S = SQ + h^A, 

where A is a self- adjoint semiclassical pseudodifferential operator of order 
zero. Since the distance between consecutive eigenvalues of 5o is 0{h), if A is 
bounded for h sufficiently small the spectrum of S forms non-overlapping 
clusters of eigenvalues 

{Ej^k; k = !,■■■ ,mj}, j = 0,l,---, 

defined by the condition 

\E,^kih)-Ejih)\ = oih^). 

As we will see the boundedness assumption of A can be dropped, provided 
one looks at the spectrum of S locally: For all A in standard symbol classes 
the spectrum of S restricted to compact intervals forms spectral clusters for 
small enough h. 

By analyzing these spectral clusters using wave-trace techniques, one gen- 
erates a collection of spectral invariants which are known as "band invari- 
ants" and have been studied in detail in various contexts in |CVj . |Weij . 
[Gulj . |Gu2| . |Dozj . |Urj . [UrV] . |Zel] . and elsewhere. The main topic of this 
paper will be computations of these invariants for the harmonic oscillator. 
These computations will be valid in arbitrary dimension, but in dimension 
two we will be able to extract from them a number of new inverse results. 
In particular, for the Schrodinger-like operator Sq + h?V, we will show that 
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if V is real analytic and symmetric in its xi and X2 coordinates, and if, 
in addition, the two eigenvalues of the Hessian of V at the origin are dis- 
tinct, then V is spectrally determined. In fact we will show more generally 
that a result of this nature is true for real analytic semi- classical potentials, 
V = V{x,h), and for C°° semi-classical potentials provided that V{x,0) is 
a quadratic form with distinct eigenvalues. (See section 6 for details.) 

Perturbations of the form Sq + h}~^^A, 6 > and A a bounded zeroth 
order pseudodifferential operator, have been studied recently by D. Ojeda- 
Valencia and C. Villegas-Blas , |Q- VV-B] . Working in Bargmann space they 
derive the first band invariant for such perturbations. 

2. The Clustering Phenomenon 

We will be working with semiclassical pseudodifferential operators, with 
amplitudes in the symbol classes and using the Weyl calculus 

(see |Marj for definitions). 



2.1. Existence of Spectral Clusters. Let 5 be as in (1.2). The following 
theorem is a special case of a result by B. Helffer and D. Robert in |HRj 
(Theorem 3.9). It states that, if we only look at eigenvalues contained in a 
fixed "window", then as ^ — )• 0, the eigenvalues of S form clusters around 
the unperturbed eigenvalues. We include a proof here for completeness, and 
to set up the computation of the band invariants. 

Theorem 2.1. For any compact interval / C M there exists ho > small 
enough such that for < h < hg, the spectrum of S in I consists of eigen- 
values that cluster near the eigenvalues of Sq within a distance 0{h?). More 
precisely, given I 

3C > 0, ^0 > such that V/i e (0, ^o] 

(2.1) Spec{S)r\ic U [hj - Ch^ , hj + Ch^] . 

Proof. Define the time-dependent operator R{t) by the identity 

(2.2) ^ih-HSo^-ih-Hs ^ J ^ 

The left hand side of this identity is a semiclassical pseudodifferential oper- 
ator (a composition of semiclassical FIOs with inverse canonical relations), 
so R{t) is a semiclassical pseudodifferential operator. In what follows we let 

R = R{2tt) = i(e-2-^'»"'^-/). 

Now let X ^ Cq^ (M) be a cut-off function which is identically equal to one 
in an open neighborhood of /, and define 

Rx = x{S)R. 
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This is a semiclassical pseudodifferential operator, and in the expansion of 
its symbol 

one can see that each rj has compact support contained in i/Q"^(supp(x)), 
where 

(2.3) H^{x,p)='^-{\x\^ + \p\'') 

is the principal symbol of Sq. It follows from the C alder on- Vaillancourt 
theorem that a bounded operator in L^, and the bound is uniform in 
h. Thus, for h sufficiently small the spectrum of / + hR^^ is contained in a 
neighborhood of 1 away from the origin in C, and we can define 

(2.4) w=:^\og[I + hR^), 

which is a bounded pseudodifferential operator of order zero. Define now 

5^ := 5 - n^W. 

Note that all the operators in this paragraph are functions of S, and therefore 
commute with each other. 

It follows from the definitions that 

g-2«r.-is^ = (/ + hR) {I + hR^)-\ 

Now let A G / be an eigenvalue of S with normalized eigenvector ip. Since 
it is clear that 

which is to say that Sy-ip = hj ip for some integer j. From this and the fact 
that S = + h'^W it follows that 

(2.5) A = hj + h^fi, where Wip = nip. 
Therefore 

\A-hj\ < Ch^, 

where C is a uniform LP' bound on the norm of W for h small enough. This 
completes the proof. □ 



Note that the intervals on the right-hand side of (2.1 ) are pairwise disjoint 
for h small enough, since their centers are a distance h apart and their lengths 
are 0{h^). The eigenvalues of S in any one of those disjoint intervals must 
be of the form 

(2.6) Ej^k = Ej + h'^fij^k 

where the fij^k are the eigenvalues of W restricted to the eigenspace of 
corresponding to Ej = Hj. Note that the fij^k are uniformly bounded. 
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Remark 2.2. The previous theorem holds for more general Zoll operators. 
Namely, let Q be a semiclassical pseudodifferential operator such that Spec(Q) C 
hZ. (This implies that the Hamilton flow of its principal symbol is 2tt pe- 
riodic.) Let A be a zeroth order semiclassical pseudodifferential operator 
whose symbol lies in 0{{x,(,)"^), and consider the perturbation of Q 

P = Q + h^A. 

Then for any compact interval / C M the eigenvalues of P in I cluster near 
hZ, for all h small enough. 



Clusters for potentials of quadratic growth. 

In the case that the perturbation A is a "multiplication by a potential 
function operator, there is a much simpler and direct proof of the 

previous result provided V is of no more than quadratic growth at infinity. 
Moreover, in this case one can see the spectral clustering appearing uni- 
formly in all of M, instead only locally uniformly in compact intervals. In 
other words, there exists C > and ho > such that for all h < ho, 

(2.7) Spec(So + h^V) C (j [hk- Ch^, hk + Ch"^] . 

Here is a sketch of the proof: If V is of no more than quadratic growth at 
infinity, then one can find constants Ci > 0, C2 > such that 

-C7i|xp -C2< V{x) < C7i|xp + C2. 

On the other hand, by directly computations, one can show that the eigen- 
values of 

So + h^i±Ci\x\'^ ± C2) = -\^^^ + ^kP(i ± Ih^Cx) ± h^C2 

are precisely hs/X it 2h?C\{^f -|- ^) it ^^C2, which are within 0{h'^) distance 
of the unperturbed eigenvalues (1.1). Now (2.7) follows from a min-max 
argument. 

Remark 2.3. Analogous arguments show that the same result holds for a 
"semiclassical potential" , 

A = Vo{x) + hVi{x) + h^V2{x) + ■■■ , 

where all the Vj are of functions with uniform quadratic growth : | Vj (x) \ < 
Ci + C2\x\'^ (for some constants Ci,C2 independent of j). 
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2.2. Cluster projectors. Let / C M be a compact interval. Keeping the 



notations of Theorem 2.1, let E € Int(/), -E > 0, and restrict the values of 
h to the sequence 

(2.8) h=^, N = l,2,.... 
For each such h small enough, 

(2.9) Spec(5) n[E -Ch^ , E + Ch^] = {E + h'^fiN,k ; k = l,... , miv} 
is a single cluster (see (|2.6[)). 



Theorem 2.4. For each h as in (2.8), let 11^ he the orthogonal projec- 
tion from L^(M") onto the span of the eigenf unctions of S corresponding 



to eigenvalues in the individual cluster (2.9). Then the family {n^j/vr is a 



semiclassical Fourier integral operator associated with the canonical relation 
{{x,y) £ M^"xM^" I Ho{x) = E = Ho{y),x lies on the same Hq orbit as y}. 

Proof. Take a closed interval Ii containing E in its interior and included in 
Int(/), and a cutoff function xi such that xi = 1 on -^i and Xi = outside 
/. Define 

U{t)=xi{S^)e-'^'^''''-. 
Then one can check that (recalling that h = E /N) , 

(2.10) ^ j^J uity'^-'^^dt = xi(5x)n^- 

Now U{t) is a Fourier integral operator associated with the graph of the 
principal symbol of . The conclusion of the theorem follows from this and 



the calculus of FIOs applied to the left-hand side of (2.10). □ 



3. Band invariants 
3.1. The first band invariant. 



According to equation (2.2), the time derivative of R{t) is given by 



(3.1) R{t) = ^[Sq, R{t)] - iA - ihAR{t). 

It follows that its principal symbol, ro(t), satisfies 

fo{t) = {Ho,ro} - ia and ro\t=o = 0, 
where a is the principal symbol of A. The solution to this problem is 

ro{t,x,p) = -i a{<ps{x,p))ds, 







where (j)s is the Hamilton flow of Hq. 
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In what follows we will denote, for any function b on M^", 
(3.2) := ^ l\{Mx,p))ds. 

With this notation, the principal symbol of 



The first band invariant arises as follows: 
Theorem 3.1. For any f e C^(M) and ip G C°°(M), the integral 

(3.3) / / (\.^t±\Pl) ^(a-) dxdp 



is a spectral invariant for the family of operators (1.2). 



Proof. Choose the interval / in Theorem 2.1 such that supp/ CC I. Then 



J2 fiEj) = TV ifiS^MW)] ■ 

j k=l 

As ^ — 7- 0, the trace on the right-hand side has an asymptotic expansion 
with leading coefficient (3.3). □ 

As a consequence, one can see that for any E > 0, the numbers 
(3.4) [ (^(a"™)dA 

J\x\^ + \p\2=E 

are spectral invariants of the semiclassical family of operators S = So + f\^A, 
where d\ is the (normalized) standard Lebesgue measure on the sphere 
|xp + = E. In fact, these numbers arise as Szego limits: 

Theorem 3.2. Fix E > and ip G C°°(M). Then, as N -f oo 



(3-5) V95(m,fc) = / 



(^(a™^)(iA + 0(l/iV). 



The proof is standard: The left hand side of (3.5) is the normalized trace 
of ip{W)Yi^. The asymptotic behavior of this trace can be computed sym- 



bolically given that the cluster projectors are FIOs (Theorem 2.4). 
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3.2. Some properties of the averaging procedure. 

For future reference we gather here some properties of the averaging pro- 
cedure, restricted to functions V of the position variable x alone. For such 
functions, letting z = x + ip, one can rewrite the integral (3.2) as 

From this integral formula one can derive many properties of V^""^ as a 
function on C". Here is a very incomplete list: 

(1) is always an even function. Moreover, one can regard y'^™ as a 
function defined on S^^'^-^/S^ = CP"-^ 

(2) corresponding complex sphere |zp = r^. As a consequence. If Vi(x) = 
V2{x) on a ball |xp < r^, then Vi^'^ = VJ™ on the complex ball 
|zp < r^. 

(3) If V{x) is an odd function, then = 0. (In later applications, we 
will mainly focus on the V^s that are even in each variable. ) 

(4) If V{x) is of polynomial growth, 

\V{x)\ < Ci + Cslxr, 



< Ci+C2\zr. 

More generally, if < ^(l^^l) for some increasing function h, 

then < h{\z\). 

(4a) In general, one cannot drop the "increasing" assumption on h 
above. For example, consider n = 1 and V{x) = jqij^, then 

V{x) = 0{j^) as X — )• oc, and 

1 r^- 1 



I/^™(z) = — / dt 



1 



27r ^ 

dt. 



2tt In 1 ^ M 



l + ^(l+cos(2t)) 

One can check V^^z) ^0{j^). 

(4b) However, if V{x) = o(l), then V{z) = o(l). We will leave the 

proof as an exercise. 
(4c) Similarly, if V{x) > C|x|^ for |x| large enough, then y^™(z) > 

C'l^l^ for \z\ large enough. 
(5) If V{x) is homogeneous of degree m in x, then V^^°{z) is also homo- 
geneous of degree m in z. 
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(6) If V{x) = x'^ = x"^ ■ ■ ■ is a monomial, then is zero unless 
|a| is even, in which case 



(3.6) 1 ^ (ai 



2i"i ^ , , Vii 

JlH hjn^^ 



^1 ^1 



(7) is identically zero if and only if V is an odd function. We will 
prove this in the n = 2 case in §5. The case n > 2 reduces to the 
two dimensional the trajectories of the harmonic oscillator 

in any dimension lie on two dimensional planes. 

3.3. Higher band invariants for perturbations by potentials. 

To derive higher order band invariants, one must look at higher order 
terms in the symbol of W. Recall that by the semiclassical Weyl calculus, 
the Weyl symbol of a composition PQ is 

a#nb ~ ef (^^^«-^^^'')(a(x,p)6(y,(?))|,=,.,,=p 

for some bi-differential operators Bj, where 

Bo(a,b) = ab, and Bi(a,b) = — {a,b}. 

2i 

In particular, the Weyl symbol of the commutator |[S'o,-R] is given by 

{Hq, r} + h^{HQ,r}^ + h^{Ho, rjs + • • • , 
where { , } is the Poisson bracket, and { , }fc is the "higher poisson bracket" 

\a\+m=k 

Let us now assume that the operator A is the operator "multiplication by 
a potential function V{x)" , i.e. 

(3.7) S = So + h^V. 

Since Hq = |(|xp + |pp) is quadratic in both x and p, all higher order 
Poisson brackets of Hq with r{t) vanish. If we write 

r{t,h,x,p) ~ ro{t,x,p) + hri{t,x,p) + h'^r2{t,x,p) H , 

and if abbreviate rk{t,x,p) to rk{t), then from ( |3.1[ ) we get 

k k j,l 
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with initial conditions rjt(O) = 0. In particular, 

ri{t) = {Ho,ri} - iVro and ri\t=o = 0. 
It is easy to check that the solution to this equation is simply 

n{t,x,p) = KQ{t,x,p). 

Similarly the equation for r2 is 

r2{t) = {Ho,r2}-i{Vn + Bi{V,ro)) and n\t=o = 0, 
and its solution is 



1 3 ^ 



6 " 2 



t rt-s 



^0 



r2it,x,p) = -r^ + - / {V{<ps{x,p)),V{<ps+u{x,p))} duds. 



In general, for any k >1, the rk{t,x,p) is the solution to the problem 

fc-i 

(3.8) rk{t) = {Ho,rk} - i'^Bi{V,rk-i-i) and rk\t=o = 0- 

1=0 

Note that if we let 

k-l 
1=0 

then gk depends only on ro, • • • ,rk-i, in particular it is independent of r^. 
So the problem (3.8) can be solved via Duhamel's principle iteratively: 

rk{t,x,p) = -i gk{t - s,4>s{x,p))ds. 
Jo 

Replacing t by 2tt, we get the higher-order symbols of R, and thus the symbol 
of Ry^ in the region Hq^{I). 

To calculate the higher-order symbols of W, note that 

w = ^^io,ii + nR,) 

So if we let 

W Wq + hwi + ■ ■ ■ 

be the full symbol of W, then wq = y*^™, and 
fc+i , 



j=l /iH h/j-i+miH \-mj=k-j+l 

where r^-'s are symbols of R^ calculated above, all indices lj,mj^s are non- 
negative, and 

S'i'-'^J-i(rmi, • • • ,rmj) = Bu_^{Bi {■ ■ ■ {Bi^{rm:,,rm2),rm3) ■ ■ ■ ),rm,_^,rm,). 
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In particular, 
(3.9) 



wi=ri- -rg = 



and 



W2 



27r 



r2 - rori + g'^o 



1 
1 

47r 



2lT r2TT-S 



{V{(ps{x,p)),V{(ps+u{x,p))} duds 



JO 

2n ru 



Jo 



{V{(f)six,p)),V{(f)uix,p))} dsdu. 



In what follows we will let, for any function F on phase space. 



(3.10) F^{x,p) :-- 

so that 
(3.11) 



1 



2tt pu 



Jo 



{F{(l)six,p)),F{(j)uix,p))} dsdu, 



W2=V' 



We are now in a position to extend Theorem 3.1 and obtain higher order 
terms in the asymptotic expansion of 

rrij 

(3.12) J2 /(^i) E = Tr [fiSMW)] • 



k=l 



Theorem 3.3. Let f G C^(M) and ip{s) = s'+^ in \3.1^j , where I > is 
an integer. Then (3.12) has a semiclassical expansion in powers of h, and 
the second coefficient is equal to 



a + /( 



I X I ^ _j_ I ?Kil 2 



^^yaveyyA ^xdp + Qf , 



f 



where Qj is equal to the integral over R^" of the following expression: 

(3.13) B2{f{Ho), (y--)'+i) + /(Fo)(y™^)W + 

+ [f'{H^){V + y''-) - nU){Ho)] . 

Here B2 is the operator in the Moyal product, we have let 

i-i 



w 



3=0 



and TZ is the operator 



nf){s) = lns) + y"\s). 



The proof is a calculation that we have sketched in the appendix. 
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3.4. The case of odd potentials. Recall that V is odd iff is identi- 
cally zero. By equations (3.9, 3.11), in this case the operator W is of order 
—2 and its principal symbol is 

Theorem 3.4. For odd potentials V, the integrals 

f{^-^^^^)^{V'^)dxdp 
where f E C^(M) and (p G C°°(M) are spectral invariants. 



Proof. Analogously as in Theorem |3.1[ the previous quantity is the coeffi- 
cient of the leading order term in the asymptotic expansion of 

rrij 

J2 fiEj) V^ih-^ l^k,j) = Tr [f{S^Mh-''W)] . 

j k=l 

In the present case h~'^W is a pseudodifferential operator of order zero and 
principal symbol V^. □ 

3.5. Perturbations by semiclassical potentials. 

More generally, we can consider the harmonic oscillator perturbed by a 
semiclassical potential, that is, consider 

(3.14) S = So + h'^V{x,h), 

where 

vix, h) ~ Vo{x) + hVi{x) + h^V2{x) + • • • 

as symbols in 0{{x, S)^). The first band variant was calculated in §3.1, with 
a{x,p) = Vo{x). For the higher invariants, a similar calculation as in §3.3 
shows that for A: > 1, 

f-fc = -{HQ,rk} -iVk-i ^ Bi{Vm,rn) and rk\t=o = 0. 

m+n+l+l=k 

Again these equations can be solved via Duhamel's principle. In particular, 
we get 

1 /■* 

n{t,x,p) = -r^ +i Vi{(t>s{x,p))ds, 

which leads to wi = V^^^. More generally, 

rk{t, x,p) = i / Vk{(t>s{x,p))ds + terms depending only on Vq, • • • , 
Jo 

and 

Wk = Vfc'™ -|- terms depending only on Vq, - ■ ■ , Vk-i. 



k-l 
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So if we take ip{x) = x'^^, then the li?^'^^ term in 

rrij 

j fe=i 

is 

f{Ho){VD^Vr + terms depending only on Fq, • • • , Vk-i. 
We conclude: 

Theorem 3.5. From the k^^ term of the expansion one can spectrally de- 
termine the quantities 

(3.15) J /( l^l' + ^P^^ ){vr)'vrdxdp + Q{{Vo, ■ ■ ■ ,Vk-i), 

where (Vq, • • • , V^-i) depends only on Vq, ■ ■ ■ , Vfe_i (and their derivatives) 
and on f, I. 

4. First inverse spectral results 

4.1. Spectral rigidity. 

Obviously for any O G SO(n), the rotated potential 

V^ix) := V{Ox) 

is isospectral with the potential V{x). From this observation one can easily 
construct trivial families of isospectral potentials. 

Definition 4.1. We say a potential V{x) is spectrally rigid if for any smooth 
family of isospectral potentials F*(x), with V'^{x) = V{x), there is a smooth 
family of orthogonal matrices Ot G SO{n), such that V*{x) = V{Otx). 



Remark 4.2. Consider the Taylor expansion of V at the origin: 

I0)x- I V 
dxj ' ^ dxidx 



V{x) ~ V^(0) + — (0)xi + — — - — (O)xiXj + higher order terms. 



It is obvious that if we rotate V{x) to V'^{x), the following remain un- 
changed: 

• the constant term V{0), 

. the length |Vy(0)|2 = ^(||(0))2, 

• the eigenvalues of the quadratic form ^ dx dx 

We will see that these data are spectrally determined in many cases. More- 
over, it is easy to see that the effect that the constant term V{0) and the 
linear term ^ §^i^)xi have on the spectrum are simply translations, and 
therefore one may drop these terms and only consider potentials of the form 

V{x) = aiX^ + higher order terms. 
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(For a semiclassical potential, one may assume the leading term Vq{x) if of 
this form. ) 

Remark 4.3. For a semiclassical potential one has another method to con- 
struct trivial isospectral families. Namely, suppose 

vix, h) = Vo{x) + hVi{x) + h^V2{x) + ■ • • 

is a semiclassical potential. Then a change of variables 



converts the operator 



Xi + hih^ 



^{-h^A + \x\^) + h^v{x,h) 



into 



2 

Therefore the semiclassical potential 

,6^ 



v{x,h) = v{x + hh^,h) + ^xibi + h^^ 

is isospectral with V{x,h). Note that if we take hi = — f^(0), then the 
leading term olV{x, h) is 



which has no linear term. Iteratively using this method, i.e. making changes 
of variables 

Xi y Xi -\- b'lX , 

m > 2, one can convert V(x,h) to an isospectral semiclassical potential 

V{x, h) = Vo{x) + hVi{x) + fev2{x) + ■■■ 

where each Vi{x) has no linear term. 

Now let V{x) be an even potential function. We say that V is formally 
spectrally rigid if for any smooth family of isospectral even potentials, V*{x), 
with V'^{x) = V{x), we have 



for all A; > 1. 



= 0. 

t=o 
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Theorem 4.4. Let V{x) be a potential which is even in each variable. Then 
V is formally spectrally rigid if it satisfies: the only even function W such 
that 



(4.2) 



/(|z|2)(^(y™^)VF™^dz(iz = 



holds for all f G C^(M) and Lp £ C°°(M) is the zero function, W = {). 

Proof. Let V^{x) be a smooth family of isospectral even potentials, such 
that F°(x) = V{x). Then 



dt 



t=o 



f{\z\^)v ((y*)^™) dzdz 

d 



f{\z\')ip'ivn 



dt 



{VT^'dzdz, 



t=o 



which implies ^\t=oV^ = 0. Similarly 

d^ 







dt^ 



t=o 



/(|z|2)(^((y*)-) dzdz 



/(l^l 



if" (y°) 



t\a,vc 



t=0 



/ /TT-avc^ 



dt^ 



t\avc 



t=0 



dzdz, 



which implies (4.1) for k = 2. Continuing in the same fashion, one can prove 
( [4T| ) for all k. □ 

Similarly one can define the formally spectral rigidity for semiclassical 
potentials 

vix, h) = Vo{x) + hVi{x) + h^V2{x) + ■■■ , 

where each Vk{x) is even, by requiring that for any smooth familhy of isospec- 
tral semiclassical potentials V^{x, h), one has 

d'' 



dt^ 



t=o 



vl = o 



for all /c > 1 and all j > 0. And one can show by using Theorem 3.4 that 
V{x,h) is formally spectral rigid if 

(1) Vq is formally spectrally rigid, e.g. it satisfies the hypotheses of 
Theorem 4.4 above. 

(2) For each j > 1, Vj satisfies the following condition: the only even 
function W such that 



/(|zn(y'=)W"™(izdz = 



holds for all / G Cg°(M) and / G N is the zero function. 
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4.2. Recovering one dimensional even potentials. 



According to (3.6), if V is an analytic function of one variable, then 
after averaging the only terms that survive are the even terms in the Taylor 
expansion of V. It follows that if n = 1 the semiclassical spectrum of 
S determines the even part of V for analytic functions V. The following 
theorem shows that the analyticity assumption can be dropped: 

Theorem 4.5. For n = 1, the first band invariant determines the even part 
ofV, i.e. V{x) + V{—x), for any smooth perturbation V . 

Proof. It is sufficient to prove the statement for V even. Then from the first 
band invariant gives us, for any r > 0, the integral 

Tv/2 

v{rcose) de. 



Making the change of variables s = r cos 6 and u = s'^, and denoting Vi{x) = 
V{y/x)/y/x, the above integral becomes 

2 2 

V?" — s Jo — u yu Jq 



The latter one equals 



r(^)(J^Vi)(r-'), 



where J 2 is the fractional derivative of order ^. So if we apply J 2 again 
and integrate, we can recover Vi, and thus V itself. □ 

Remark 4.6. In |GUWj . we showed that by using semiclassical invariants 
modulo 0(^^), one can spectrally determine not only V{x) + V{—x), but 
also V'^{x) + V'^{—x). It follows that can determine V{x) itself under suit- 
able symmetry conditions. It turns out the same result holds in higher 
dimensions, i.e. one can spectrally determine the integral / V{x)dax as well 
as the integral f V^{x)dax over any sphere \x\ = r. In particular, one can 
distinguish radially symmetric potentials from other potentials. For more 
details, c.f. |GUW) . 

4.3. Recovering one dimensional odd analytic potentials. Now as- 
sume y is a one dimensional odd analytic potential, 

V{x) = aix + a^x^ 

Then 
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From the first band invariant for odd potentials (see Theorem 3.4), 



we get, by replacing fi hy ^ and making change of variables z — )• Az, the 

spectrally determined expression 

(4.3) 



/• /'27r ru 

7~^j Jc Jo Jo 



k,l odd 

Note that in complex coordinates, 

fr . c-dfdg 
^^^'^ = ^d-zd-r 

so we have for = 1 and / = 2Zi + 1 odd, 

f>27r 



I g-MkP I " r{e''z + e-''z, (e^"z + e-^"z)'}dsdn dzdz 
Jc Jo Jo 



Mkl^ r Qle'^'-^'^e'-^z + e-'^zy-^dsdudzdz 



^ |2 



C JO JO 

/-I 



^ ^ 27r /'li 

/( ' ] I e'^'^^''^'^\z\''^^dzdz I / sm{s-u)dsdu 

[l - l)/2J Jq Jq Jq 

1\ hi 



Similarly for k = 2ki + 1 and / = 1 we have 

io Jo A;i!|i^i+^ 



In particular, we get from the lowest order term of A in (4.3) (i.e. the coef- 
ficient of A^) the number ^y-af from which one can recover ai up to a sign. 
At this point we must distinguish two cases. 

Case 1: ai ^ 0. Suppose we have recovered the first m coefficients 
d-i,''' )fl2m-i) up to an overall sign (the same for all). Then by looking 



at the coefficient of A^™"*"^ in the expression (4.3) one can recover the num- 
ber aia2m+i + fl2m+ifli, siuce all other products of pairs are known precisely. 
It follows that a2m+i is determined up to the same sign ambiguity as oi. 

Case 2: ai = 0. More generally, suppose we have already found that 
ai = ■ ■ ■ = ak-2 = 0, where /c > 3 is an odd number. Let I > k he odd too. 
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Then a calculation as above shows that 

/ e-^l^l' r{{e''z + e-''z)^,{e'''z + e-'''zy}dsdudzdz 
Jc Jo Jo 

^ k - 2m\k - 1 - mj y-=J^ J' 

m=0 ^ / \ 2 / 



where C is some non-vanishing constant that depends on /j,. Note that 

1 /A: - 1\ - 1' 
k\k-l, 



l-k 
2 



> 



and for < m < |, 

so the expression 

/ e-'^l^l' r{{e''z + e-''z)\ (e*"z + e-'^'zyjdsdu dzdz 
Jc Jo Jo 

is non- vanishing as well. It follows that for I = k, one can determine a^, and 
for Z > A;, one can determine a^ai, from which one can determine ai up to 
the same sign ambiguity as ak- 
in conclusion, we have proved: 

Theorem 4.7. For any odd analytic potential V , in one dimension, the first 
band invariant for odd potentials determines V{x) up to a sign. 

4.4. Recovering constant and quadratic terms. 

In all dimensions V{0) is spectrally determined, since 

F(0) = V^^iO) = hm / y^™dA. 

'■^0j|a;|2+|p|2=r2 

About recovering the quadratic terms, we have already seen that wc can 
only hope to recover the eigenvalues of the Hessian at the origin. Since the 
average of any odd function vanishes, we can assume that V{x) is even itself. 
We will assume moreover that V{x) is analytic, so that up to a rotation, 

V{x) = V{0) + ^aixj + V{x), 

where V{x) contains only terms that are homogeneous in x of degree greater 
than four. It follows that 

= F(0) + J2 (^Mf + V^'^'iz), 
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where V^^^^z) contains only terms that are homogeneous in z of degree 4 
and higher. So from the spectral invariant 

J e-^'^'^'vJ^+Uzdz 

we get, by replacing /x by ^ and making change of variables Zi — >■ Xzi, the 
expressions 

// \ k-\-l 

e-MNP (\2 ai\zin dzdz + terms involves at least A2'=+^ 

It follows that the integral 

is spectrally determined for all k. In particular, one can determine 

[ e-'^N'eS»'l^'l'dzdz = 7r"^ 

J jJL-ai ii-an 

and thus determine the polynomial 

(/X - ai) - an). 

This proves: 

Theorem 4.8. For analytic potentials the eigenvalues of the Hessian at the 
origin are spectrally determined. 

4.5. Recovering the linear term for certain potentials. First assume 
V is an analytic odd potential, so that 

where Vk is homogeneous of degree k, and in particular, V\{x) = W{Q) ■ x. 
Then according to Theorem 3.4, the integrals 

are spectrally determined. Again we replace /x by ^ and make change of 
variables z — )■ Az, wc get the expression 

J e~''l^l dzdz + terms involves at least 

So one can spectrally determine the integral 

from which one can read off ||Vy(0)|p. 
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One can improve this result slightly by considering analytic potentials of 
the form 

(4.4) V{x) = Vq{x) + Vi{x) + y3(^) + y^{x) + Vq{x) + --- , 

i.e. analytic potentials with vanishing quadratic and quartic terms. Note 
that this property, 

^2 = and = 0, 

is a spectrally property: We have already seen in section 4.4 that one can 
determine whether V2 = using spectral data. In this class of potentials, 
one can spectrally determine the integrals 

{vr)d\. 



ISr 

If all these integrals vanish, then we get 



^itv-dX = Yo\{Sr), 



avc 



from which one can deduce that V?"" = 0, and thus V4 = 



Now suppose V{x) is a potential of the form (4.4). Without loss of gen- 
erality, one can assume that the constant term is zero. Then according to 
theorem 3.3, with / = and f{t) = e~^^, one can spectrally determine 



dzdz. 



Since V2^^ = V^**™ = 0, after replacing // by and making change of 
variables 2: — t- Az, the lowest order term is again 

and from this one gets ||Vy(0)|p. In conclusion, we have 

Theorem 4.9. Let V{x) be an analytic potential. Then one can spectrally 
determine whether V{x) is of the form |^.^| ), and ifV is in this class, one 
can spectrally determine its linear term Vi(x) up to rotation. 

5. The first band invariant in dimension 2 

The averaging procedure on the space of functions 

C^^^iR^) := {Vixi,X2) = V{xlxl), V G C^iR^)} 

of smooth functions on which are even in each variable has particularly 
nice properties that we investigate in this section. 

Let TT be the cotangent fibration vr : = T*R? R?, and let be 
the average of tt*V with respect to the circle action 

(5.1) e'^z = ie'^zi,e'^Z2). 
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As we have seen, for ip E C°°(M), the integrals 

(5.2) Vi.iy)= /'e-2(l^ll' + l^2p)^(yave)^^^- 



are spectral invariants of the operator (1.2). 

To analyze these invariants we'll begin by decomposing into its 

Fourier coefficients with respect to the circle action 

(5.3) e''z = {e''zi,e"'z2), 

and this we'll do by first examining the one-dimensioanl analogues of these 
Fourier coefficients. More explicitly, let 

(5.4) Br : ,(M) ^ n(K) 
be the operator 



f{x)=f\^ — J =g(|z|V-)^5r, 

where gr is the 2r-th Fourier coefficent of g{\z\'^e^'^^) with respect to 6. By 
definition this operator commutes with the homothety x ^ Ax, A G M+, and 
hence maps x^^ into a multiple, 'yk,rx'^^, of itself. To compute 7^^^ we note 
that 

-\ 2fc r=k / r)T \ 

Z -\- Z \ 1 I \ , ,2k„2ire 



r=—k 



Hence 



for |r| < k and zero for |r| > k. Note that Bq is invertible. 

To describe the asymptotic dependence of 'yk,r on for A; >> we apply 
Stirling's formula 



n 



I ~ V27m (-) (l + ain ^ + a2n ) 



to the quotient 

If 2fc \ _ {2k)\ 
A''\k + rJ {k + r)\{k-ry. 



V^47r2(fc2 - r2)(Mr)fc+4(fc^)fc-r- ^ ^ 

-L=(l + ci,^A;-i + •••), 
V irk 
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giving the asymptotic expansion 



^ oo 

(5.6) 7fc,r ~ y^c;,,.A;"', co,r = 1, 
or (more intrinsically) the asymptotic expansion 

oo 

(5.7) Br = ^Ck,rQ''-'^ 



r=0 



where Q = f Hence in particular is a pseudodifferential operator on 
Ce^en(IK) of Order 



2- 

Coming back to the problem in two dimensions that prompted these com- 

lOO / 

even \ 



putations, for each V G C^g„(M^) let RnV the nth Fourier coefficient of V 



with respect to the circle action ( |5.3[ ). We will prove: 
Theorem 5.1. The operator 

(5.8) Rn : n(M') ^ C'e™n(M') 

is zero for n ^ mod 4, anrf modulo the identification C^^^ = C^^^{M.)'S'C^^j^ 

(5.9) Rn = Br^Br, u = Ar 

(i.e. for functions of the form /(xi, X2) = g{xi)h{x2),Rnf = Bj.g{xi)Brh{x2) ). 
In particular Rq is invertible. 

Proof. It suffices to check this for g = xf^ and h = x"^ in which case 
yave ^ JL4-(fe+0 f (^^gie ^ zie-^^)2^(z2e^^ + Z2e-'^f^ dO 



AA'AtJ^'^' ^^^^ 

|r|<min(fc,i) 

where, in polar coordinates zj = |zj|e*^^ , 

(5.10) Vk,i,r = 4-('=+') (^/_^^^) |zi|2'=|z2p^e2-(^-^^). 

Hence setting Xj = \zj\ we get for the n-th Fourier coefficient, n = 4r, of 
ynvc ^[^i^ respect to the circle action (5.3), 

hk,rx\)i%rxi) = BrX^Brxl 

□ 



e 

2k\ /2r 
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Remark 5.2. The operator Ar = Rmn = 4r, can also be described as a 
Radon transform. Namely consider the double fibration 



M X M M X M 

where 7r( 21,^2) = {Rezi,Rez2) and p{zi, Z2) = (|zi|, |2;2|)- Then : C^cni^'^) 
C^gn(M^) is the transform 

(5.11) ArV = P*(^j rl^e,^*Ve-^'''' dOidO^ , 

where Te^^eM, Z2) = e^^i(e^^2^i, e-^^2^2). 



If we take ip{t) = in (5.1) we get the spectral invariants 

(5.12) j e-2(l^ll' + l^2p)(yavc)m^^^- 

and if we expand y^^° in its Fourier series with respect to the circle action 



(5.3), i.e. express as the sum 

ArV {xi,X2)e 



Arid 



then (5.12) becomes 



(5.13) /'e-i(-?+-i) ( Ar,V ■ ■ ■ Ar^v] dxidx2. 

In particular for m = 1 we obtain 

(5.14) J e-2(='i+''2)AoVdxidx2, 



and for m = 2 



(5.15) J e-2('-'i+'4)Y^\ArV\^dxidx2. 



6. Inverse spectral results in two dimensions 

Throughout this section we work in two dimensions, and consider the 
inverse spectral problem for a perturbation of the harmonic oscillator by a 
potential (at times semicassical). 



6.1. Results for smooth perturbations. 
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6.1.1. Potentials of the form + f2{x\). In this subsection we consider 

potentials of the form 

(6.1) V{xuX2) = fi{xl) + f2{xl). 
It is easy to see that 

(6.2) ^"™(x,p) = M^I+pD + V2{xI+pI), 
where 

r(_i/2) ,y2[fj{sy 



(see the proof of Theorem 4.5). 



Theorem 6.1. Generically, potentials of the form (6.1) are spectrally de- 
termined, up to the obvious symmetries: Exchanging of the roles of xi and 
X2, and adding a constant to fi and subtracting the same constant from f2- 
(The genericity condition is M^). 



Proof. As we have seen (c.f. (3.5)), the integrals 



(6.3) / {v''^'Y{d)de, 

where Sr C is the sphere of radius r, are spectrahy determined. We claim 
that the previous integral equals 

(6.4) 27rV / [ipi[r'^{l-u)) + ^2{r'^u)\'' du. 

Jo 

To see this, consider the map $ : — )• given by 

$(x,p) = {xj +pI, xl +pI) =: {r'^ui,r'^U2). 

This maps the sphere Sr onto the line segment ui + U2 = 1 (in the first 
quadrant). We will take u = U2 as a coordinate in that segment. The 
fibers of $ are tori, which we can parametrize in complex coordinates by 
(r-y/iiie*'^, ry^e**). The variables (s,t, u) parametrize 5^. Using these vari- 



ables to compute (6.3) yields (6.4) (we'll omit the details of the calculation). 
The integrals (6.4) determine the distribution function of 

ipr{u) = 93i(r^(l - u)) + (p2{r^u) 

on [0, 1]. We now make the genericity assumption that 

(6.5) '>Pr{u) is monotone for each r > 0. 

and therefore its distribution function determines it up to the ambiguity 
^r{^ ~ u), which amounts to switching the roles of xi and X2. Finally, it 
is not hard to see that knowing the two-variable function ip determines ipi 
and ip2 up to the ambiguity of adding and subtracting a constant. □ 
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6.1.2. Semiclassical Potentials with quadratic Vq. 
Theorem 6.2. Semiclassical potentials of the form 

V{x, h) = ax\ + hxl + hVi + h^V2 H , 

where a^h and Vi G C^^^{M?) for each i, are spectrally determined. 



Proof. According to Theorem 3.5, one can spectrally determine 



plus a term not depending on Vi with i >1. We can rewrite this invariant 
as 

From this we can determine 

for all /z and A > sufficiently large with respect to /i (so that the integral 
is absolutely convergent). Note that the phase equals 

— (A + afi)xi — (A + h^)x\, 

and since a / 6 we can make a linear change of variables, (A, //) — >■ (/ii, ^2) 
so that the invariant becomes 

g-Mix-2-^2xi AQVidxidx2. 

Integrating in polar coordinates we get, up to a universal constant, the 
spectral invariant 



I' 



-fiir 



'-''''^'AoViiV^^,V^2)/iV^iV^2)dridr2 



This function of fii and ^2 is known spectrally in some infinite wedge of 
(determined by the condition that A > 0). However, it is easy to see 
that this is an analytic function of {111,112), and therefore it is spectrally 
determined in the whole quadrant M^. Using the inverse Laplace transform, 
one can determine AqVi{xi, X2) pointwise, and thus determine Vi itself. 

By an inductive argument, one can similarly spectrally determine each of 
the remaining Vfc. □ 



6.2. Inverse Spectral Results for Real Analytic Perturbations. 
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6.2.1. Spectral rigidity of the quadratic potentials ax\ + hx"^. In this section 
we will prove the following result, which will be used later: 

Proposition 6.3. If a ^ b the potential V = ax\ + hx\ is ^'formally" spec- 
trally rigid; i.e. if Vj, —e < t < e, is a smooth family of potentials each of 
which is an even function of xi and X2 and has the same spectrum as V , 
then 

d^ 



(6.6) 
for all k. 



dt^" 



Vt 







Proof. We first note that if ( 6.6 ) is non-zero for some k then we can reparametrize 
t so that (6.6) is non-zero for k = 1. Thus it suffices to prove that the func- 
tion W 

0, we get 



^|t=o is zero. 



To see this we first note that by remark 5.2, 



for i 7^ and = V . Thus by inserting Vt into (5.13), differen- 
tiating with respect to t and setting t 



for all m. Thus 



and since 



AqW dxidx2 = 



(xl + xl) 



-V 



a + ua 2 



x\ 



2' ' 2 2 ' 2 

and a ^ b, we can, by making a linear change of coordinates, u) 
(/Xi,/i2), convert this equation into the form 



. iii^2_M2 „2 



'^2AoW{xi,X2)dxidx2 = 



for all Hi, fJ-2 in an infinite wedge of the 1st quadrant. Arguing just as in the 



end of the proof of Theorem 6.2 we can apply the inverse Laplace transform 
and conclude that AqW = 0, and by the injectivity of Aq, that W = 0. □ 

6.2.2. Spectral determinacy. Since one can spectrally determine the integrals 
of X] l^j^P over the circles -|- x| = t, in particular one can determine 
whether this sum vanishes to infinite order at t = xi = X2 = 0, and hence 
whether V itself vanishes to infinite order at t = 0. This can be strengthened, 
as follows: 

Theorem 6.4. // V is an analytic function on such that its Taylor 
expansion at the origin is of the form 



V = V2 + Va + V6 + 
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with each Vj homogeneous of degree j, and the quadratic term V2 has two 
distinct eigenvalues a b, then V is spectrally determined (up to a rota- 
tion)^ 

Proof. We can assume without loss of generality that V2 = axf + 6x2, and 
let 

y'^™ = a|zi|2 + 6|z2p + Vr + + ■■■ 
be the Taylor series expansion of y^™. From the spectral invariants 

we get, by replacing by ^ and making change of variables Zi — )■ Xzi, the 
expressions 



^2k+4 J ^-H(|,,|2 + |,,|2) („|^^|2^5|^2|2)'=+1 

+ {k + 1)A2'=+6 y g-H(|,,|2 + |,,|2) ^^|^^|2 ^ ^|^^|2)fc y..e ^ 0{X^^+^) 

and hence, from the proof of the spectral rigidity of axf + and linearity 
of this expression in V^'™, this suffices to determine V^*™. 

Continuing: the next unknown term in the expansion above is 

ik + l)X^'+^ J e-t(l-il'+NI') {a\zi\^ + b\z2\^)''vrdzdz 

and hence, by the same argument, this determines Vg'™ (and by induction 
and repetition of this argument Vg'™, Viq^, ■ ■ ■)■ Thus the Taylor series of 



-J/ ave jg spectrally determinable and hence by (3.6) so is the Taylor series of 



V. □ 

6.2.3. The case of semiclassical potentials. Now consider perturbations of 
by the semiclassical potentials, 

+ + h\Vo + hVi + h^V2 + --- ), 

where 

(6.7) Vq = axf + + higher order terms 

with a 7^ 6 as above, and all Vi's are even functions. As we have seen from 
the previous section that the first spectral invariant is enough to determine 



^Results of this nature for the Schrodinger operator, — ft'^Au" + V,V real analytic, can 
be found in |GU] , |GPU| and [Hez| . However, these results require strong "non-rationality" 
assumptions on the coefficients of the leading term ^ aiX^ of V. 
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Vq. To determine Vi's with i > 1, we use the higher invariants. We have 
seen in section 3.5 that the integral 

is spectrally determined. If we write 

v, = vi + v^ + vi + --- 

with each V-^ homogeneous of degree i, then the same argument as in the 
previous subsection shows that 

l+Tn=k 

are spectrally determined. Since each is already known (as component of 
Vq), the above functions suffice to determine each of V^, and thus determine 
Vi. By an inductive argument, one can determine all the 1^'s. 
To summarize, we have proved: 

Theorem 6.5. A semiclassical potential of the form X]fc>o f^^Vk where Vq 



is of the form (6.1) with a ^h, and, for each k is analytic and satisfies 



Vfc(±xi, ibx2) = Vk{xi,X2), is spectrally determined. 

Appendix A. Proof of theorem 13.31 



To prove Theorem 3.3 we must compute the first two non-trivial terms in 
the full symbol of the operator f{S-^)W''~^^. 

A.l. The symbol of f{Sy^). Let as^ ~ Hq + h'^S2 + • • • be the asymptotic 
expansion of the symbol of S^. Note that S2 = V + y^™. The operator 

U{t) = e**^^ 

is a pseudodifferential operator of order zero. Let 

CTt/ ~ uo + hul + h^U2 H 

be the expansion of its symbol. Writing down the equation —iJJ = S^U 
symbolically it is easy to check that uq = e^^^°, and one computes that 

-iui = Bi{Hq,uo) + HqUi, ui\t=o = 0, 

which has as solution ni = 0. U2 is the solution to the problem 

(A.l) - iu2 = B2{Ho,uo) + H0U2 + U0S2, n2|i=o = 0. 

Using the general formula 

\a\+m=2 
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one finds: 



(A.2) B2{Ho, no) = J {t^Ho - int) e^*^". 



Substituting into (A.l) and integrating one obtains: 



(+3 
Ho + — +t{V + V^ 
12 01 

We can now prove the following 

Lemma A.l. Let (yf[Sx) ~ 4>o + h4>i + f\^ (^2 + ' ■ ■ the asymptotic expansion 
of the full symbol of f{S^). Then one has: (/)o = f{HQ), 0i = and 

11 1 

(A.4) <j>2 = {V + y-) /'(Fo) - -f"m - -Ho f"'{Ho). 

Proof. By the Fourier inversion formula, 

1 



/(5x) = ^ J U{t)f{t)dt. 

At the symbolic level this reads: 't'j = ^ f Uj{t) f{t)dt. The stated formulae 
follow from this and the previous calculations of the Uj, j = 0,1,2. □ 

A.2. The end of the proof. First we need to compute the second non- 
trivial term in the expansion of the symbol of VF^"'"^: 

Lemma A.2. For any non-negative integer I and smooth functions wq, W2 
on M^", one has: 

I 

{wo + h^W2r^'+^^ = ^0+' +h^{l + l)wiw2 + Y,wiB2{wo,W^'] 

The proof is by induction on I. 



Using lemmas A.l and A.2 and recalling that 

wo = and W2 = V^, 

one can easily compute the second non-trivial term in the full symbol of 



f{S)W^^^. The result is precisely the function (3.13) plus 



This finishes the proof of Theorem 3.3 



30 



V. GUILLEMIN, A. URIBE, AND Z. WANG 



References 

[CV] Y. Colin de Verdiere, "Sur le spectre des operateurs elliptiques a bicharateristiques 
toutes periodiques" , Comment. Math. Helv. 54 (1979), 508-522. 

[Doz] S. Dozias, "Clustering for the Spectrum of /i-pseudodifferential Operators with Pe- 
riodic Flow on an Energy Surface", J. Funct. Anal. 145 (1997), 296-311. 

[Gul] V. Guillemin, "Band asymptotics in two dimensions". Adv. Math. 42 (1981), 248- 
282. 

[Gu2] V. Guillemin, "Spectral Theory on S^: Some Open Questions", Adv. Math. 42 
(1981), 283-298. 

[GPU] V. Gillemin, T. Paul and A. Uribe. "'Bottom of the well' semi-classical trace 

invariants". Math. Res. Lett. 14 (2007), no. 4, 711-719. 
[GU] V. Guillemin and A. Uribe, "Some Inverse Spectral Results for Semi-classical 

Schrodinger Operators", Math. Res. Lett. 14 (2007), 623-632. 
[GUW] V. Guillemin, A. Uribe and Z. Wang, "A semiclassical heat trace expansion for 

the perturbed harmonic oscillator", preprint. arXiv: 1107.2960 
[GuW] V. Guillemin and Z. Wang, "Semiclassical Spectral Invariants for Schrodinger 

Operators", submitted. 
[HR] B. Helffer and D. Robert. "Puits de potentiel generalises et asymptotique semi- 

classique", Ann. Inst. Henri Poincare 41 (1984), 1-43. 
[Hez] H. Hezari. "Inverse Spectral Problems for Schrodinger Operators", Comm. Math. 

Phys. 288 (2009), 1061-1088. 
[Mar] A. Martinez, An Introduction to Semiclasscial and Microlocal Analysis, Berlin- 
Heidelberg- New York: Springer, 2002. 
[O-VV-B] D. Ojeda- Valencia and C. VlUegas Bias, "On limiting eigenvalue theorems in 

semiclassical analysis", preprint, 2011. 
[Ur] A. Uribe, "Band Invariants and Closed Trajectories on S""', Adv. Math. 58 (1985), 

285-299. 

[UrV] A. Uribe and C. Villegas-Blas, "Asymptotics of Spectral Clusters for a Perturbation 
of Hydrogen Atom", Comm. Math. Phys. 280 (2008), 123-144. 

[Wei] A. Weinstein, "Asymptotics of eigenvalue clusters for the Laplacian plus a poten- 
tial", Duke Math. J. 44 (1977), 883-892. 

[Zel] S. Zelditch, "Fine Structures of ZoU Spectra", J. Funct. Anal. 143 (1997), 415-460. 

Department of Mathematics, Massachusetts Institute of Technology, Cam- 
bridge, MA 02139, USA 

E-mail address: vwg@math.mit.edu 

Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, 
USA 

E-mail address: uribe@uinich.edu 

Department of Mathematics, University of Michigan, Ann Arbor, MI, 48109, 
USA 

E-mail address: wangzuoq@umich.edu 



